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fSJ . The Shubnikov - de Haas effect in quasi-two-dimensional normal metals is studied. The interlayer 

» I ' conductivity is calculated using the Kubo formula. The electron scattering on short-range is consid- 

CL(, ered in the self-consistent Born approximation. The result obtained differs from that derived from 

■ the Boltzmann transport equation. This difference is shown to be a general feature of conductivity 

in magnetic field. A detailed description of the two new qualitative effects - the field-dependent 
phase shift of beats and of the slow oscillations of conductivity is provided. The results obtained are 
applicable to strongly anisotropic organic metals and to other quasi-two-dimensional compounds. 
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Magnetic quantum oscillations were discovered long ago and were frequently used as a powerful tool of studying the 
geometry of Fermi surfaces and other electronic properties of various metals. In recent years, quasi-two-dimensional 
(quasi-2D) organic metals ||^ attract great interest because many new unconventional effects are very pronounced 
in these compounds. These effects are high-Tc superconductivity, spin and charge density waves, strong anisotropic 
magnetic quantum oscillations etc. Much work was devoted to studying magnetic quantum oscillations in these 
compounds (for a review see e.g. ||^). The quantum oscillations of magnetization are a thermodynamic effect that is 
completely determined by the density-of-states distribution. Any exact calculation of the electron density of states 
] (DoS) is a very complicated problem but a semi-phenomenological description of magnetization oscillations in quasi- 
' 2D compounds was recently provided in a number of theoretical papers The chemical potential oscillations and 

O [ the arbitrary electron reservoir due to the open sheets of the Fermi surface create no principal difficulties . Since 
, the number of occupied LLs is very large {np > 100) in most quasi-2D organic metals, the effect of the electron- 

electron interaction is reduced (as in the Fermi liquid) and can be taken into account via the renormalization of the 
electron effective parameters. Differences in Landau level shape (which depends on a particular compound) lead to 
^ ■ only limited quantitative differences in magnetization curves yielding no qualitatively new effects. On the qualitative 
' level, therefore, the de Haas-van Alphen (dHvA) effect in quasi-2D normal metals is believed to be well understood. 
Attempts of theoretical description of the quasi-2D Shubnikov-de Haas (SdH) effect were not as successful although 
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^ _ some work on this subject appeared in recent years |^,^-10|. There are still many open qualitative questions. 

. One of these open questions is the origin of the phase shift in the beats of the resistivity oscillations with respect to 
fSj ' those in the magnetization. The beat behavior of the oscillations in quasi-2D metals is known reliably to originate from 
a slight warping of their Fermi surfaces in the direction normal to the 2D plane. The superposition of the contributions 
from the maximum and minimum cyclotron orbits leads to an amplitude modulation of the k-th harmonic by the 
^ I factor cos(27rfcAF/2i3 — 7r/4), where B is the magnetic field and AF = {ch/2Tre){Amax — ^mm) is the difference 
' between the oscillation frequencies caused by the extreme orbits with the fc-space areas Amax and Amin , respectively 
101 . From the beat frequency one can readily evaluate the warping of the Fermi surface and hence the interlayer 
transfer integral At « ep/S.F/F (see e.g. |l^Jll[]). The situation becomes less clear when the warping is so weak that 
less than one half of the beat period can be observed experimentally. In principle, an observation of one single node 
Q , would already be quite informative ||l^, provided the phase offset (i.e. the phase of the beat at 1/ B 0) is known. 
1| In the standard Lifshitz-Kosevich (L-K) theory |13| this phase offset is strictly determined by geometrical reasons 
, and is equal to —tt/A for both the dHvA and SdH effects (l). However recent experiments on layered organic metals 
^ ■ K-(BEDT-TTF)2 Cu[N(CN)2]Br Q] and (BEDT-TTF)4[Ni(dto)2] |l| revealed a significant difference in the node 
H \ positions of the beats of dHvA and SdH signals. The phase shift in the latter compound was estimated to be as big 
. P.\ as 7r/2. 

Another very interesting phenomenon (also not explained in the framework of the standard theory) is slow oscil- 
lations of magnetoresistance that were observed in a number of quasi-2D organic metals [pHl^-[l9|] . The behavior of 
slow oscillations resembles that of the SdH effect that lead to a suggestion of additional, very small Fermi surface 
pockets in these materials. However, band structure calculations (basically giving a good description of the electron 
band structure and of the Fermi surface topology in organic metals) show no evidence of such small pockets in any of 
these compounds. Moreover, while slow oscillations are often very pronounced in magnetoresistance, sometimes even 
dominating the oscillation spectrum, thus far no analogous observation in oscillating magnetization (dHvA effect) was 
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reported. Certainly, slow oscillations carry useful information about the compounds but one needs some theoretical 
explanation and, desirably, a quantitative description of the phenomenon to extract this information. Note that in 
both cases the oscillation spectrum was strongly dominated by the first harmonic when no substantial deviations from 
the standard L-K theory are expected. 

An explanation and the qualitative description of phase shift of beats (which is based on the Boltzmann transport 
equation) was proposed recently |^ together with a comparative experimental study of this effect. An idea that 
the slow oscillation may arise as an entanglement of different rapidly oscillating contributing factors in conductivity 
which have slowly oscillating amplitudes due to beats was also recently suggested in and proved by presenting 
experimental results on temperature and angular dependences of slow oscillations. In this paper we give a more 
accurate theoretical description of these phenomena as well as a detailed calculation of the interlayer conductivity in 
quasi-2D metals in strong magnetic field. 

In sec. II the general formula (|o|) is derived. The calculation is performed starting from the Kubo formula. 
The result of this calculation differs from that obtained using the Boltzmann transport equation. The origin of this 
difference is pointed out. The additional term is shown to be a general feature of conductivity in magnetic field. 
However it becomes essential only in the quasi-2D case. The simple explicit formula for interlayer conductivity is 
obtained in the self-consistent Born approximation in sec. III. A discussion of the reliability and possible application 
of the results obtained is given in sec. IV. 



II. GENERAL FORMULA FOR INTERLAYER CONDUCTIVITY 

We consider a quasi-2D metal in magnetic field perpendicular to the conducting layers: B\\z. The electron spectrum 
of quasi-2D electron gas in magnetic field is then given by 

e{n,k^) = hujc {n + 1/2) - 2tcos{k^d) (1) 

where t is the interlayer transfer integral, kz is the wavevector perpendicular to the layers, d is the interlayer distance, 
uJc = eB/in*c is the cyclotron frequency. Both huic and t are assumed to be much smaller than the Fermi energy. 

To calculate conductivity we use the Kubo formula |^2|. The procedure is similar to that in three-dimensional 
metals without magnetic field ( |22[ , § 7.1.2). In magnetic field only the new set of quantum numbers m = {n, kz, ky} 
should be used instead of momentum p and the alternative dispersion relation (^. The evaluation of the Kubo formula 
without vertex corrections gives 



(2) 



where the volume V normalizes the sum over quantum numbers m , e is the electron charge, the limits of the integral 
over e are (— c»;(xi), n'p{e) is the derivative of the Fermi distribution function: 

- 7i'p{e) = l/{4Tcosh2 [(e - ^i)/2T]} (3) 

and A{m, e) is the spectral function that is related to the electron Green's function G^{m, e) or to the retarded 
self-energy part I]^(m, e) : 

„a, . — 2ImS^(m, e) , , 

A(m, e = -2ImG^ (m, e = ^ (4) 

[e - e(m) - Rel]«(m, e)]^ -f [ImS-R(m, e)]^ 

Formula (||) is close to the corresponding formula without magnetic field [ |2|], formula (7.1.10)] until the self-energy 
part I]^(m, e) is specified. It arises mainly from impurity scattering. The main contribution to resistivity comes from 
the short-range impurity scattering. We approximate short-range impurities by point-like ones. Then, if one does not 
take into account the diagrams with intersections of impurity lines in the self-energy (the contribution of such diagrams 
at finite kz dispersion of electrons is usually small) the electron self-energy depends only on electron energy and not 
on electron quantum numbers. This fact greatly simplifies the calculations because the sum over quantum numbers m 
in formula (^) can be now computed analytically. The constant part of the real part ReS-'^(e) of electron self-energy 
produces only a constant shift of the chemical potential. It does not influence the physical effects and, hence, is 
omitted in the subsequent calculations. The small oscillating part ReE^ (e) of ReE^(e) enters the final expression for 
conductivity in the second order in small damping factors. Hence, it can affect slow oscillations of conductivity and 
should be kept in accurate quantitative analysis. It always come in the combination e* = e — ReS^(e). 
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On the contrary, the imaginary part of self-energy ImE^(e) is very important since it describes the momentum 
relaxation of electrons. 

Performing the summation over ky in (|^) and changing integration over by integration over energy e{n,kz) we 
get 

„ 7o J 27r 



where Nj^j^ = B/^qcI is the electron density on one Landau level and the electron velocity u^(e,n) is given by 
To go further we have to transform the sum over LL s to a sum over harmonics. This can be done using the Poisson 
summation formula (Appendix A). Substituting (AS) into (^ we obtain: 



^ (-1)^ exp — Ji —— / —A^e', e) (-n'^(e)) 



27r ^^-^ hk ^ \ hujc ) ^ \ hujc ) 

°° , ,.2td'^ f ATikt\ f de ,,,,,, , 



-^AT ^ ' .^k'^td^ , fATrkt\ f de 
e Nll } ^ 



where one should use the expansion Ji{kx)/k = x/2 for the zeroth harmonic fc = 0, and the integral Iz{e,k) over e' 
can be easily evaluated with the spectral function (Q): 

-2ImI]^(e) \ f2nike' 




exp 



f2mke*\ ( 1 2-Kk\ ^ , , . 



where e* = e — ReS-'^(e) and 



[ImS«(e)]^y V 
1 27rfc 
\ hiOc ) \|ImS^(e)| fiLOc 

i?^(fc,e) = exp(-27r|fc| |lmE^(e) I /Tiw^) (9) 



has the form similar to that of the usual Dingle factor RD{k) = exp (^—2TT'^k kBTo/hLUc) ■ Collecting formulas (^ and 
(i) we get 



k— — oo 



{-lf2td'^ ^ /ATTkt\ 

\ hujc J 



'2mke*\ ( 1 27rfc\ , 

><exp^— ™?7:^ + ^ i?D(fc,e). (10) 



hujc ) \ |ImI]-^(e)| TiLJc) 

Note that this expression has additional term 2'Kk/hijJc near the standard 1/ |lmE^(e)| term in round brackets in 
the second line. This term can not be obtained from the Boltzmann transport equation (compare, for example, with 
the results of j2^] and |2^) and it arises only due to quantization of electron energy spectrum in magnetic field. This 
quantization results in fast oscillations of the mean square electron velocity as function of energy which give rise to 
the rapidly oscillating factor exp {2'Kike' / fiuj c) in (^). This factor appears only in magnetic field. The derivative of 
this factor with respect to e' comes after the integration over e' in (0) because the function A^{e',e) — {G^ — G^Y 
has one second-order pole in each complex half-plane. These second-order poles arise from the combinations [G^Y 
and (G^)^, where G^ and G^ are the advanced and retarded Green's functions respectively. So, the additional term 
in (nOI) has quantum origin; it appears due to fast oscillations of the mean squared velocity at the Fermi level in 
magnetic field. 

To go farther, we need an explicit form of the electron self-energy which enters formula (p^). It is calculated is in 
the next section. 
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III. CONDUCTIVITY IN SELF-CONSISTENT BORN APPROXIMATION 



We consider electron scattering only by short-range impurities because these impurities make the main contribution 
to the relaxation of electron momentum. To calculate the electron self-energy we use the self-consistent Born ap- 
proximation. The graphical representation of the Dyson equation for the irreducible self-energy part in self-consistent 
Born approximation is shown in fig. ^ By such approximation we neglect multiple scattering on one impurity (no 
more than two dash lines go to one impurity in fig. The single dash line in fig. |l| corresponds to the first-order 
term which leads only to a constant shift of the chemical potential and, hence, can be omitted. 

a a 

f A 

I / \ 

FIG. 1. The Dyson equation for the irreducible self-energy in self-consistent Born approximation. The double solid line 
symbolizes exact electron Green's function. 

The corresponding analytical expression is 

^"{m,e) ^ l^U^ G{n,n,E)^ ^ C,U^ j d^rG{r,r,E) (11) 

where is a sum over all impurities and the brackets (..) denote averaging over impurity positions, Ci is concentration 
of impurities which are assumed to be uniformly distributed p^ . The electron Green's function G'(r, r, E) in formula 
( ITll ) contains the self-energy determined by the same formula (^l|) (this is why the approximation (|l|) is called 
self-consistent Born approximation). The Green's function is isotropic along the conducting planes. Hence, one can 
write 

G{r,r,E) = \c^{z)\'G{E), (12) 

where the electron wave function 4i{z) along z-axis does not enter the final result because it disappears after integration 
over z in (fo), and 



G{E) = i A{E) + in 



l + 2g(-l)^Jo(^jexp(^2...^-^ 



(13) 



A{E) is a slowly varying function of energy which can be taken at the Fermi energy. The exact form of this function 
is not important for conductivity in Born approximation. Formula ( p^ ) can be derived performing the summation 
over electron quantum numbers m = {n, kz,kx} in the definition of the Green's function: 

where En^k^,k^ is the electron energy in the state with quantum numbers m = {n, kz,kx} given by (^. The electron 
wave function ^'n.fcj,fc^(^) in Landau gauge is approximately given by 

*n,fc,,fe,(r) = — j=^^xniy - yo)ct>{z) 

where i/q = —cfikx/eB and the normalization condition |Xn(2/)|^ dy = 1 is used to perform integration over kx- 
The further calculation of the sum in (14) is similar to that in (||). 



The Born approximation (formula (^T|)) takes into account only the first term of expansion in the small parameter 
ttUNll/Tt-^c = T^f Id, where Nll/TilOc is equal to the electron density of states at the Fermi level in unit volume, / is 
the scattering amplitude (which is constant at small wave vector q <C 1/ro, Tq is the range of the impurity potential). 
For short-range impurities the parameter / /d is usually small. 

From (|l^) one can easily see that in Born approximation the imaginary part of self-energy is proportional to the 
density of states pa : 
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- ImS^(e) = ~aU^ ImG'(e) = tt dU^ x p{e). (15) 

The unknown coefficient C = tt CiU^ in (|l5|) is simply related to the average Dingle temperature Tq: (|lmE^(TO, e)|) = 
C ■ (/3(e)) — C ■ {Nll/TiuJc) (1 + Tin) — TrfesT^), where the triangular brackets mean an average value of a quantity 
inside, fc^ — 1-38 • IQ^^^erg / K is the Boltzmann constant and nn is the density of reservoir states that exist in many 
organic metals due to the open sheets of the FS. 

In the extreme 2D case {Tuoc 3> t), substantial deviations from formula (|5|) are possible because the strong 
degeneracy of the LLs makes the Born approximation not applicable. Since we consider the case 2t > hujc (when 
the beats of the oscillations exist) and f /d <^ 1, we shall use (^ij) for our subsequent calculations that are now 
straightforward. 

From the formulas (|ll], |l| and [l|) we have 

|lml]^(m,e)| = irksTD ( 1 + 2f;(-l)'=Jo { cos { i?i3(fc,e) ) . (16) 

V fe=i 



Together with (^) this gives a nonlinear equation for ImS^ (rn,e). We can solve it in the strong harmonic damping 
limit using iteration procedure i.e. by making an expansion in the small oscillating part which is an expansion in a 

parameter (^RD\/fiuJc/2Tr'^t^ . To treat the slow oscillation accurately one has also to pick up all second-order slowly 



oscillating terms. We get 



|ImS^*(e)| « TrkBTD-{ 1 - 2 Jq ( ^] cos ] Rod } ■ (17) 



where Rqd — exp {~2Tr'^ kBTo/huJc). There is no slowly oscillating second-order term in the self energy in Born ap- 
proximation. At this point the real part of the electron self energy was important because it canceled the contribution 
from the entanglement with the oscillations of the Dingle factor (|9[). In the second order in damping factors any 
combination of the form 

f 2Tr{e-Ret^{e))\ / -27r |lml]«(e) I \ /27re\ „ 

cos exp '— = cos - — Rod 

y nuJc J y nujc J \ nujc J 

does not produce slowly oscillating term. This statement can be easily checked by substituting ( pT| ) with (^3|) into 
(|l6|). If we neglected ReE^(e) in (|l^) we would get an additional slowly oscillating term 
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in curly brackets of (17) which arise from the mixing with the oscillating Dingle factor and enters not only the 
imaginary part of the self-energy but also the density of electron states (see (p^). The slow oscillations of p{e) 
would result in huge slow oscillations of magnetization which are increased by an additional factor ep/htOc- Such 
huge slow oscillations of magnetization have surely not been observed and the result (|l^) is correct. Beyond the 
Born approximation the relation (|l^) is no more valid and Iml]^(e) acquires some slow oscillating term although 
p(e) does not. In the previous (unpublished) work |2^ ] the real part of electron self-energy was disregarded and 
the entanglement with the oscillations of Dingle temperature (producing additional large contribution to the slow 
oscillations) was incorrectly taken into account. 

Substituting (|l^ into (|l^) we obtain the following expression for the conductivity: 



1 - 



2^o(|S:)cos(lf)i?z,(e) 



t ^ \ TuOr 



cos 



/27re* 



Roie) 
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If the transfer integral is large enough, Ant > Tu^c, one can use the expansions of the Bessel function at large value 
of argument: 



Jq{x) w \j2/nx cos {x — 7r/4) , x ^ 1 
Ji{x) w ^J2/ttx sin [x — tt / 4) , x ^ 1 . 



(18) 



Then performing again expansion in small parameter (^R jy ^ huJc/ ^ir'^tj and making use of the standard trigonometric 
formulas we get 



Ozz = -T -. { 1 + 2a/ — cos I I X 



fi uJcnkBTo 

f Ant TT , 
X cos I — - + 4>b] RdRt 



2ttH 



Rn* \/l + cos 



Ant n 



(19) 



where the phase shift of beats is 



(pi, = arctan (a) ; a = - — 1 H 

2nt \ rkOr 



and the phase of slow oscillations is 

05 = arctan (05) /2 where 05 = hujc/2nt. 



(20) 



(21) 



The star index in the Dingle factor Rd- of slow oscillations is explained in the next section. The temperature smearing 
factor is given by the usual L-K expression: 



Rt = 



2n'^kBT/hu 



sinh {2n^kBT/huJc) ' 



It appear in the fast Shubnikov oscillations after integration over energy of a rapidly oscillating function of energy with 
the Fermi distribution function. The slowly oscillating term depends only on the transfer integral t and is independent 
of energy. Hence, it does not acquire any temperature smearing. 

The phase shift (20) obtained from the Kubo formula is larger than that of [ formula (9)] obtained using 



the Boltzmann transport equation by a factor (l + 2n'^kBTu/huJc) ■ This difference comes from the additional term 
2nk/huc near 1/ |lmE^(e)| in round brackets in the second line of ([lo|). As has been noted after formula (|o|), this 
term has quantum origin. However, the result ( [T9| ) and ( pTf ) concerning the slow oscillation does not differ from [ pl| , 
formula (4)]. Thus the Boltzmann transport equation and the Kubo formula in self-consistent Born approximation 
give the same amplitude and phase of slow oscillations. 



IV. DISCUSSION OF THE RESULTS 



In this paper a detailed calculation of the interlayer magnetotransport in quasi-2D normal metals is performed. 
The specific features of quasi-two-dimensionality and strong magnetic field result in pronounced qualitative effects 
such as the phase shift of beats of conductivity oscillations and the slow oscillations which cannot be described 
in the framework of the standard three-dimensional theory usually applied to quasi-2D compounds. The beats of 
magnetoresistance oscillations in layered compounds are used for estimating the interlayer transfer integral which 
strongly influences different electronic properties of strongly anisotropic compounds. The field-dependent phase shift 
of beats may lead to the errors in this estimate. Hence, a detailed quantitative description of this phenomenon is 
important. 

The result of the calculation using the Kubo formula is different from that obtained using Boltzmann transport 
equation. An additional term in conductivity (see formula (|lO|) and the discussion after this formula) is general for 
conductivity in magnetic field. It arises due to Landau quantization of electron energy spectrum. However, this term 
is proportional to the oscillations of the electron mean square velocity which are smaller than the oscillations of the 
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electron relaxation time (or of — ImS^(e) ^ p(e)) by the same factor hujc/^-Kt as the phase shift of beats is. Therefore, 
this additional term in 3D case is smaller than the main oscillating term by the factor of TiiiJc/f-F- 

The slow oscillations of magnetoresistance can give useful information about compounds under study. A signifi- 
cant feature of the slow oscillations is that their Dingle factor is different from the factor Re, of the Shubnikov 
oscillations. The usual Dingle factor includes all temperature-independent mechanisms of smearing of fast quantum 
oscillations. These are not only microscopical scattering events of electrons but also macroscopic spatial inhomo- 
geneities of the sample. These inhomogeneities lead to macroscopic spatial variations of the electron energy e* in 
formula ([l0| ) which is equivalent to a local shift of the chemical potential. The total signal is an average over the 
entire sample and such macroscopic inhomogeneities lead to the damping of magnetic quantum oscillations similar 
to that caused by temperature. Since the slow oscillations do not depend on /i, they are not affected by this type 
of smearing and the corresponding Dingle temperature of slow oscillations is determined by only short-range 
scatterers. One can therefore estimate relative contributions from macroscopic inhomogeneities and from local defects 
to the scattering rate by comparing Td and T£,. This role could be quite essential in organic metals. For example, 
such a comparison for a sample of /3-(BEDT-TTF)2lBr2 gives ^ (0.8 ± 0.02) K while T^, = (0.15 ± 0.02) K, 

that means that the long-range crystal imperfections are important for damping of fast quantum oscillations. The 
relaxation of electron momentum affecting transport quantities comes, however, mainly from short-range impurities. 

The slow oscillations in formula ( |l9|) do not have temperature damping factor. Hence, although the amplitude of 
the slow oscillations contains the square of the Dingle factor (they are a second-order effect), it can be larger than 
the amplitude of the fast SdH oscillations at T > Tq. However, our experience disagrees with the statement that 
the slow oscillations do not manifest any temperature dependence (and, hence, could be seen at room temperature). 
Actually they do have some temperature damping because the oscillating DoS itself has some temperature dependence. 
The temperature damping of the DoS oscillations comes from the electron-phonon and electron-electron interactions. 
In normal 3D metals the electron-electron (e-e) scattering rate l/ree ~ {kBT)^ /hfi while the electron-phonon 
scattering rate 1/Tph ^ {kBT/K) {kBT/hiOn)^ ■ One can estimate the effect of these scattering processes on the DoS 
oscillations by introducing the additional damping factor 

RtD ~ exp [-7r(l/WcTee + l/UcTph)] (22) 

analogous to the usual Dingle factor. This factor enters squared in the amplitude of slow oscillations. The temperature 
Tso at which the slow oscillations become damped by this factor is much higher than the characteristic temperature 
of the damping of fast quantum oscillations. It is approximately given by 7r(l/a;cTee('7s'o) -I- 1 / tiJcTphiTso)) ~ 1- The 
above analysis of the temperature dependence of slow oscillations is very approximate. A rigorous calculation must 
be based on the exact calculation of the electron self-energy due to these two types of interactions. Nevertheless, the 
above arguments can give qualitative estimates. A more accurate calculation would be useful since the temperature 
dependence of slow oscillations at high enough temperature may give additional information about the electron-phonon 
and electron-electron interactions in various compounds where slow oscillations exist. This is important since these 
interactions in layered organic metals determine the superconducting and the density-wave transitions. 

The entanglement with the oscillations of chemical potential contributes an additional tem perature-dependent term 



to the slow oscillations of conductivity. This term can be easily obtained by substituting (B2) into (|19D. However, this 
term has additional damping factors and (i?£)/^D*)^compared to the main slowly oscillating term. Therefore, 
this correction is as small as the second harmonic of Shubnikov oscillations is, and we can neglect it. 

The slow oscillations does not appear in magnetization because there is no suitable entanglement of different 
oscillating quantities in magnetization. The magnetization being a thermodynamic quantity is completely determined 
by the electron density of states. However, the density of states does not have slowly oscillating terms. The mixing 
with the oscillations of the chemical potential, or with those of the Dingle factor and of ReE^(e) does not also lead 
to slow oscillations of magnetization (see Appendix B). 

Now we shall discuss the approximations made during the derivation of formula (p^). 

The first limitation of the proposed analysis is that the magnetic field is taken to be perpendicular to the conducting 
layers. A finite tilt angle 9 of the magnetic field with respect to the normal to the conducting planes may be 
approximately taken into account by rescaling the Landau level separation, lUc — > tOc cos 9, and of the warping of 
the Fermi surface [^, t{9) = t(0) Jo(fci?dtan0), where kp is the in-plane Fermi momentum. But this is only a 
semiclassical approximation based on the assumption that the FS remains the same. Actually, the tilting of the 
magnetic field changes the dispersion relation and a more profound study of the effect of tilting of magnetic field on 
transport properties is required. The quantum mechanical calculation of the dispersion relation in tilted magnetic 
field in the first order of the transfer integral gives |^ t{9)/t{0) — exp (— g^/4) L° (g^/2), where g = di&ii9 / an , 
an = ^/fic/eBz is the magnetic length and (x) is Laguerre polynomial. This result is also approximate, but it 
should work satisfactory at not too great tilt angles. In the limit n — > oo the two above results coincide. 
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More essential errors may come from the so-called incoherent or weakly incoherent electron interlayer transport. 
This means, that when the interlayer transfer integral becomes comparable to the Dingle temperature, the electron 
interlayer jumps an electron scattering on impurities should not be considered separately. Their entanglement may 
result even in qualitatively new phenomena in the quasi-2D magnetotransport In the present work we consider 
only the case t > Td. 

Other errors may come from the approximate expression for self-energy (|l^). The self-consistent Born approximation 
at finite kz dispersion works quite well, but other scattering mechanisms (especially for the calculation of the DoS) 
should be taken into account. An accurate study of this problem may depend on a particular type of the compound 
in hand. It may lead to some quantitative modifications of formula (|l9|). 

The above analysis does not take into account the vertex corrections. In our case (of point-like impurity scattering) 
this is right because, according to the Ward identity, the vertex r(m, E) — p + ni Vpl]^(m, E). Hence, if the retarded 
self-energy depends only on energy, the vertex corrections are zero. The fact that E^(m, e) is approximately a function 
of energy e only is a consequence of the short-range (or point-like) impurity potential. More precisely, if one takes a 
point-like impurity potential and neglects all diagrams with the intersections of the impurity lines in the self-energy, 
then after averaging over randomly and uniformly distributed impurity positions one obtains S^(m, e) = E^(e). 
The neglected graphs with the intersections of the impurity lines describe the coherent scattering on two impurities 
simultaneously. The contribution of such scattering is small at large enough interlayer transfer integral. In the 
three-dimensional case without magnetic field the vertex corrections produce an additional factor (1 — cos a) in the 
integrand for the transport scattering relaxation time where a is the scattering angle. But the scattering probability 
is independent of the scattering angle in the case of point-like impurities and the additive cos a vanishes after the 
integration over angles. Hence, the vertex corrections vanish. 

In derivation of formula (^^ only first- and second-order terms in the small damping factors Rt and Rd were 
taken into account, assuming the harmonic damping to be strong. This is valid in the most experiments on quasi-2D 
organic metals where the amplitude of the second harmonic does not usually exceed 5% of the first harmonic. 

So, in spite of the approximations made in the above analysis, the proposed theoretical description is valid in a 
large domain of parameters which one has in real experiments on quasi-2D organic metals. For example, above result 
concerning the slow oscillations and the phase of beats of the SdH oscillations is important for an experimental study of 
these effects in /3-(BEDT-TTF)2lBr2 (see, e.g. |0). The proposed results may be also used for heterostructures 
with large enough interlayer jumping. 

The author thanks M.V. Kartsovnik, W. Biberacher, A.M. Dyugaev and I. Vagner for encouragement and stimulat- 
ing discussions. The work was supported by the EU ICN contract HPRI-CT-1999-40013 and RFBR No. 00-02-17729a. 

APPENDIX A: TRANSFORMATION OF THE SUMS OVER LLS TO THE SUMS OVER HARMONICS 

To transform the sums over LL number into the harmonic sums we shall apply the Poisson summation formula p^ ] 

s2"'="/(n)dn (Al) 



where a G (no — 1; nn). This formula is valid for arbitrary function f{n). The electron velocity is determined from the 
dispersion relation (Q) as 

de{n,kz) 2td . 

Vz{e,n) = — — = —sm{kzd} = 

nokz n 



n 



- {t-nujc{n + ll2)f. (A2) 



The sum in (^) now becomes 



n 71=0 V \ \ 

7 °° POO 
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^^-^ h k ^ \ huJc ) ^ \ flLOc ) 

In this formula for the zeroth harmonic fc = one should use the expansion Ji[kx)/k = x/2. 



APPENDIX B: MAGNETIZATION 

The first harmonic of the oscillating part of magnetization is given by (see , formula 6) 

where the oscillating part of the chemical potential is ( Q, formula 5) 

7r(l + nfl(e_F)) V TvjJc J \nujc J 

The entanglement of magnetization oscillations with the oscillations of the Dingle factor produces an additional 
term 



(X sin — Jo - — X cos — Mo \ t 

\ riLUc ) \ TlLOc J \ If^c J \ iT^c 



which give rise to the second harmonic but makes zero contribution to the slow oscillations of magnetization 



The entanglement with the oscillations of the chemical potential (B2) produces the term 

'27 



oc sm 



3in ('^^\ 



f2'Kil{B)\ 



/'2neF\ /27r£_F\ / 47r A 2RtRd 



which also contribute only to the second harmonics (or higher harmonics) but not to slow oscillations. 



[1] Shoenberg D. "Magnetic oscillations in metals", Cambridge University Press 1984 

[2] T. Ishiguro, K. Yamaji and G. Saito, Organic Superconductors, 2nd Edition, Springer- Verlag, Berlin, 1998. 

[3] J. Wosnitza, Fermi Surfaces of Low- Dimensional Organic Metals and Superconductors (Springer- Verlag, Berlin, 1996); 

M.V. Kartsovnik and V.N. Laukhin, J. Phys. I France 6, 1753 (1996); J. Singleton, Rep. Prog. Phys. 63, 1111 (2000). 
[4] V.M. Gvozdikov, Fiz. Tverd. Tela 26, p. 2574 (1984) [Sov. Phys. Solid State 26 (9), p. 1560 (1984)]; K. Jauregui, V.I. 

Marchenko, I.D. Vagner, Phys. Rev. B 41„ 12922 (1990); M.A. Itskovsky, T. Maniv and I.D. Vagner, Phys. Rev. B 61, 

14616 (2000); P. Grigoriev and I. Vagner, Pis'ma Zh. Eksp. Teor. Fiz., 69, 139 (1999) [JETP Letters 69, 156 (1999)]; T. 

Champel and V.P. Mineev, Phil. Magazine B 81, 55 (2001); T. Champel, Phys. Rev. B 64, 54407 (2001). 
[5] N. Harrison et al., Phys. Rev. B 54, 9977 (1996). 

[6] P. Grigoriev, Zh. Eksp. Teor. Fiz. 119(6), 1257 (2001) [JETP 92, 1090 (2001)]. 
[7] A.E.Datars and J.E.Sipe, Phys. Rev. B 51, 4312 (1995) 



9 



[8] P. Moses and R.H. McKenzie, Phys. Rev. B 60, 7998 (1999) 

[9] V.M. Gvozdikov, Fiz. Nizk. Temp. 27 (9/10), p. 956 (2001) [Sov. J. Low Temp. Phys. 27(9/10) (2001)]. 
[10] N.S. Averkiev, L.E. Golub, S.A. Tarasenko and M. Willander, J. of Phys. Cond. Matter 13, 2517 (2001). 
[11] A.P. Mackenzie, S.R. Julian, A.J. Diver et al., Phys. Rev. Lett. 76, 3786 (1996); C. Bergemann, S.R. Julian, A.P. Mackenzie, 

S. NishiZaki and Y. Maeno, Phys. Rev. Lett. 84, 2662 (2000). 
[12] J. Wosnitza, G.W. Crabtree, K.D. Carlson, H.H. Wang and J.M. Williams, Physica B 194-196, 2007 (1994); J. Wosnitza, 

G. GoU, D. Beckmann, S. Wanka, D. Schweitzer and W. Strunz, J. Phys. I France 6, 1597 (1996). 
[13] A. A. Abrikosov, Fundamentals of the theory of metals, North-Holland, 1988 

[14] H. Weiss, M. V. Kartsovnik, W. Biberacher, E. Bahhes, A.G.M. Jansen, and N.D. Kushch, Phys. Rev. B 60, R16259 
(1999). 

[15] M.Schiller, W.Schmidt, E.Bahhes et al., Europhys. Lett. 51, 82 (2000). 

[16] M. Kartsovnik, V. Laukhin, V. Nizhankovskii and A. Ignat'ev, Pis'ma Zh. Eksp. Teor. Fiz. 47, 302 (1988) [Sov. Phys. 

JETP Lett. 47, 363 (1988)]; M. Kartsovnik, P. Kononovich, V. Laukhin and I. Schegolev, Pis'ma Zh. Eksp. Teor. Fiz. 

48, 498 (1988) [Sov. Phys. JETP Lett. 48, 541 (1988)]; M. Kartsovnik, P. Kononovich, V. Laukhin, S. Pesotskii and 

LF. Schegolev, Pis'ma Zh. Eksp. Teor. Fiz. 49, 453 (1989) [Sov. Phys. JETP Lett. 49, 519 (1989)]; M. Kartsovnik, V. 

Laukhin and S. Pesotskii, Fiz. Nizk. Temp. 18, 22 (1992) [Sov. J. Low Temp. Phys. 18 13 (1992)]. 
[17] T. Terashima, S. Uji, H. Aoki, M. Tamura, M. Kinoshita and M. Tokumoto, Solid State Commun. 91, 595 (1994); E. 

Ohmichi, H. Ito, T. Ishiguro, G. Saito and T. Komatsu, Phys. Rev. B 57, 7481 (1998). 
[18] B. Narymbetov, N. Kushch, L. Zorina et al., Eur. Phys. J. B 5, 179 (1998); T. Togonidze, M. Kartsovnik, J. Perenboom, 

N. Kushch and H. Kobayashi, Physica B 294-295 435 (2001). 
[19] L. Balicas, J. Brooks, K. Storr, D . Graf, S. Uji, H. Shinagawa, E. Ojima, H. Fujiwara, H. Kobayashi, A. Kobayashi and 

M. Tokumoto, |cond-mat/0008287| (2000). 
[20] P.D. Grigoriev, M.V. Kartsovnik, W. Biberacher, N.D. Kushch, P. Wyder, Phys. Rev. B 65, 60403(R) (2002) 



[21] M.V. Kartsov nik, P.D. Grigoriev, W. Biberacher, N.D. Kushch, P. Wyder, submitted to Phys. Rev. Lett, (cond- 



mat/0203279) (2002) 



[22] G. Mahan "Many-Particle Physics", 2nd ed.. Plenum Press, New York, 1990, (§7.1) 
[23] "Standard Mathematical Tables and Formulae", CRC Press, (1996) 

[24] In real layered compounds the impurity distribution may have periodic modulation in the direction perpendicular to the 

layers and with period equal to the interlayer distance. However, in Born approximation this leads only to a constant factor 

in (^l|) which renormalizes the Dingle temperature. 
[25] This fact has evident physical meaning. The matrix element of scattering by a point-like impurity U is the same for all 

transitions not changing electron energy. Hence, the scattering rate 1/T(e) — — 2ImE^/?!, is proportional to the number of 

states to which an electron may scatter, that is, to the density o f states at a given e nergy. 
[26] P.D. Grigoriev, M.V. Kartsovnik, W. Biberacher and P. Wyder, |cMd-mat/0108"352| (2001) (unpublished). 
[27] This fact can be understood in a different way: slow oscillations originate from local electron motion and are not sensitive 

to global variations of the position of the chemical potential with respect to the bottom of the conducting band. 
[28] Yamaji K., J. Phys. Soc. Jpn. 58, 1520 (1989). 
[29] Yasunari Kurihara, J. Phys. Soc. Jpn. 61, 975 (1992). 



10 



